Accurate, automatic determination of astigmatism and phase with Ctfplotter in IMOD.

Ctfplotter in the IMOD software package is a flexible program for determination of CTF parameters in tilt series images. It uses a novel approach to find astigmatism by measuring defocus in one-dimensional power spectra rotationally averaged over a series of restricted angular ranges. Comparisons with Ctffind, Gctf, and Warp show that Ctfplotter's estimated astigmatism is generally more reliable than that found by these programs that fit CTF parameters to two-dimensional power spectra, especially at higher tilt angles. In addition to that intrinsic advantage, Ctfplotter can reduce the variability in astigmatism estimates further by summing results over multiple tilt angles (typically 5), while still finding defocus for each individual image. Its fitting strategy also produces better phase estimates. The program now includes features for tuning the sampling of the power spectrum so that it is well-represented for analysis, and for determining an appropriate fitting range that can vary with tilt angle. It can thus be used automatically in a variety of situations, not just for fitting tilt series, and has been integrated into the SerialEM acquisition software for real-time determination of focus and astigmatism.

Software tools for automated transmission electron microscopy.

The demand for high-throughput data collection in electron microscopy is increasing for applications in structural and cellular biology. Here we present a combination of software tools that enable automated acquisition guided by image analysis for a variety of transmission electron microscopy acquisition schemes. SerialEM controls microscopes and detectors and can trigger automated tasks at multiple positions with high flexibility. Py-EM interfaces with SerialEM to enact specimen-specific image-analysis pipelines that enable feedback microscopy. As example applications, we demonstrate dose reduction in cryo-electron microscopy experiments, fully automated acquisition of every cell in a plastic section and automated targeting on serial sections for 3D volume imaging across multiple grids.

Building bridges between cellular and molecular structural biology.

The integration of cellular and molecular structural data is key to understanding the function of macromolecular assemblies and complexes in their in vivo context. Here we report on the outcomes of a workshop that discussed how to integrate structural data from a range of public archives. The workshop identified two main priorities: the development of tools and file formats to support segmentation (that is, the decomposition of a three-dimensional volume into regions that can be associated with defined objects), and the development of tools to support the annotation of biological structures.

Automated tilt series alignment and tomographic reconstruction in IMOD.

Automated tomographic reconstruction is now possible in the IMOD software package, including the merging of tomograms taken around two orthogonal axes. Several developments enable the production of high-quality tomograms. When using fiducial markers for alignment, the markers to be tracked through the series are chosen automatically; if there is an excess of markers available, a well-distributed subset is selected that is most likely to track well. Marker positions are refined by applying an edge-enhancing Sobel filter, which results in a 20% improvement in alignment error for plastic-embedded samples and 10% for frozen-hydrated samples. Robust fitting, in which outlying points are given less or no weight in computing the fitting error, is used to obtain an alignment solution, so that aberrant points from the automated tracking can have little effect on the alignment. When merging two dual-axis tomograms, the alignment between them is refined from correlations between local patches; a measure of structure was developed so that patches with insufficient structure to give accurate correlations can now be excluded automatically. We have also developed a script for running all steps in the reconstruction process with a flexible mechanism for setting parameters, and we have added a user interface for batch processing of tilt series to the Etomo program in IMOD. Batch processing is fully compatible with interactive processing and can increase efficiency even when the automation is not fully successful, because users can focus their effort on the steps that require manual intervention.

Practical Experience with Hole-Free Phase Plates for Cryo Electron Microscopy.

Phase plate (PP) imaging has proven to be valuable in transmission cryo electron microscopy of unstained, native-state biological specimens. Many PP types have been described, however until the recent implementation of the "hole-free" phase plate (HFPP), imaging has been challenging. We found the HFPP to be simple to construct and to set up in the transmission electron microscopy, but care in implementing automated data collection is needed. Performance may be variable, both initially and over time, thus it is important to monitor and evaluate image quality by observing the power spectrum. We found that while some HFPPs gave transfer to high resolution without CTF oscillation, most reached high resolution when operated with modest defocus.

MRC2014: Extensions to the MRC format header for electron cryo-microscopy and tomography.

The MRC binary file format is widely used in the three-dimensional electron microscopy field for storing image and volume data. Files contain a header which describes the kind of data held, together with other important metadata. In response to advances in electron microscopy techniques, a number of variants to the file format have emerged which contain useful additional data, but which limit interoperability between different software packages. Following extensive discussions, the authors, who represent leading software packages in the field, propose a set of extensions to the MRC format standard designed to accommodate these variants, while restoring interoperability. The MRC format is equivalent to the map format used in the CCP4 suite for macromolecular crystallography, and the proposal also maintains interoperability with crystallography software. This Technical Note describes the proposed extensions, and serves as a reference for the standard.

The challenge of determining handedness in electron tomography and the use of DNA origami gold nanoparticle helices as molecular standards.

The apparent handedness of an EM-tomography reconstruction depends on a number of conventions and can be confused in many ways. As the number of different hardware and software combinations being used for electron tomography continue to climb, and the reconstructions being produced reach higher and higher resolutions, the need to verify the hand of the results has increased. Here we enumerate various steps in a typical tomography experiment that affect handedness and show that DNA origami gold nanoparticle helices can be used as convenient and fail-safe handedness standards.

Clustering and variance maps for cryo-electron tomography using wedge-masked differences.

Cryo-electron tomography provides 3D imaging of frozen hydrated biological samples with nanometer resolution. Reconstructed volumes suffer from low signal-to-noise-ratio (SNR)(1) and artifacts caused by systematically missing tomographic data. Both problems can be overcome by combining multiple subvolumes with varying orientations, assuming they contain identical structures. Clustering (unsupervised classification) is required to ensure or verify population homogeneity, but this process is complicated by the problems of poor SNR and missing data, the factors that led to consideration of multiple subvolumes in the first place. Here, we describe a new approach to clustering and variance mapping in the face of these difficulties. The combined subvolume is taken as an estimate of the true subvolume, and the effect of missing data is computed for individual subvolumes. Clustering and variance mapping then proceed based on differences between expected and observed subvolumes. We show that this new method is faster and more accurate than two current, widely used techniques.

Exploring the retinal connectome.

PURPOSE: A connectome is a comprehensive description of synaptic connectivity for a neural domain. Our goal was to produce a connectome data set for the inner plexiform layer of the mammalian retina. This paper describes our first retinal connectome, validates the method, and provides key initial findings. METHODS: We acquired and assembled a 16.5 terabyte connectome data set RC1 for the rabbit retina at ≈ 2 nm resolution using automated transmission electron microscope imaging, automated mosaicking, and automated volume registration. RC1 represents a column of tissue 0.25 mm in diameter, spanning the inner nuclear, inner plexiform, and ganglion cell layers. To enhance ultrastructural tracing, we included molecular markers for 4-aminobutyrate (GABA), glutamate, glycine, taurine, glutamine, and the in vivo activity marker, 1-amino-4-guanidobutane. This enabled us to distinguish GABAergic and glycinergic amacrine cells; to identify ON bipolar cells coupled to glycinergic cells; and to discriminate different kinds of bipolar, amacrine, and ganglion cells based on their molecular signatures and activity. The data set was explored and annotated with Viking, our multiuser navigation tool. Annotations were exported to additional applications to render cells, visualize network graphs, and query the database. RESULTS: Exploration of RC1 showed that the 2 nm resolution readily recapitulated well known connections and revealed several new features of retinal organization: (1) The well known AII amacrine cell pathway displayed more complexity than previously reported, with no less than 17 distinct signaling modes, including ribbon synapse inputs from OFF bipolar cells, wide-field ON cone bipolar cells and rod bipolar cells, and extensive input from cone-pathway amacrine cells. (2) The axons of most cone bipolar cells formed a distinct signal integration compartment, with ON cone bipolar cell axonal synapses targeting diverse cell types. Both ON and OFF bipolar cells receive axonal veto synapses. (3) Chains of conventional synapses were very common, with intercalated glycinergic-GABAergic chains and very long chains associated with starburst amacrine cells. Glycinergic amacrine cells clearly play a major role in ON-OFF crossover inhibition. (4) Molecular and excitation mapping clearly segregates ultrastructurally defined bipolar cell groups into different response clusters. (5) Finally, low-resolution electron or optical imaging cannot reliably map synaptic connections by process geometry, as adjacency without synaptic contact is abundant in the retina. Only direct visualization of synapses and gap junctions suffices. CONCLUSIONS: Connectome assembly and analysis using conventional transmission electron microscopy is now practical for network discovery. Our surveys of volume RC1 demonstrate that previously studied systems such as the AII amacrine cell network involve more network motifs than previously known. The AII network, primarily considered a scotopic pathway, clearly derives ribbon synapse input from photopic ON and OFF cone bipolar cell networks and extensive photopic GABAergic amacrine cell inputs. Further, bipolar cells show extensive inputs and outputs along their axons, similar to multistratified nonmammalian bipolar cells. Physiologic evidence of significant ON-OFF channel crossover is strongly supported by our anatomic data, showing alternating glycine-to-GABA paths. Long chains of amacrine cell networks likely arise from homocellular GABAergic synapses between starburst amacrine cells. Deeper analysis of RC1 offers the opportunity for more complete descriptions of specific networks.

Electron tomography of RabA4b- and PI-4Kß1-labeled trans Golgi network compartments in Arabidopsis.

The trans Golgi network (TGN) of plant cells sorts and packages Golgi products into secretory (SV) and clathrin-coated (CCV) vesicles. We have analyzed of TGN cisternae in Arabidopsis root meristem cells by cell fractionation and electron microscopy/tomography to establish reliable criteria for identifying TGN cisternae in plant cells, and to define their functional attributes. Transformation of a trans Golgi cisterna into a Golgi-associated TGN cisterna begins with cisternal peeling, the formation of SV buds outside the plane of the cisterna and a 30-35% reduction in cisternal membrane area. Free TGN compartments are defined as cisternae that have detached from the Golgi to become independent organelles. Golgi-associated and free TGN compartments, but not trans Golgi cisternae, bind anti-RabA4b and anti-phosphatidylinositol-4 kinase (PI-4K) antibodies. RabA4b and PI-4Kß1 localize to budding SVs in the TGN and to SVs en route to the cell surface. SV and CCV release occurs simultaneously via cisternal fragmentation, which typically yields ∼30 vesicles and one to four residual cisternal fragments. Early endosomal markers, VHA-a1-green fluorescent protein (GFP) and SYP61-cyan fluorescent protein (CFP), colocalized with RabA4b in TGN cisternae, suggesting that the secretory and endocytic pathways converge at the TGN. pi4k1/pi4k2 knockout mutant plants produce SVs with highly variable sizes indicating that PI-4Kß1/2 regulates SV size.

ER sliding dynamics and ER-mitochondrial contacts occur on acetylated microtubules.

The endoplasmic reticulum (ER) network is extremely dynamic in animal cells, yet little is known about the mechanism and function of its movements. The most common ER dynamic, termed ER sliding, involves ER tubule extension along stable microtubules (MTs). In this study, we show that ER sliding occurs on nocodazole-resistant MTs that are posttranslationally modified by acetylation. We demonstrate that high MT curvature is a good indicator of MT acetylation and show in live cells that ER sliding occurs predominantly on these curved, acetylated MTs. Furthermore, increasing MT acetylation by drug treatment increases the frequency of ER sliding. One purpose of the ER sliding on modified MT tracts could be to regulate its interorganelle contacts. We find that all mitochondria and many endosomes maintain contact with the ER despite the movements of each. However, mitochondria, but not endosomes, preferentially localize to acetylated MTs. Thus, different ER dynamics may occur on distinct MT populations to establish or maintain contacts with different organelles.

Cryo-electron tomography of microtubule-kinesin motor complexes.

Microtubules complexed with molecular motors of the kinesin family or non-motor microtubule associated proteins (MAPs) such as tau or EB1 have been the subject of cryo-electron microcopy based 3-D studies for several years. Most of these studies that targeted complexes with intact microtubules have been carried out by helical 3-D reconstruction, while few were analyzed by single particle approaches or from 2-D crystalline arrays. Helical reconstruction of microtubule-MAP or motor complexes has been extremely successful but by definition, all helical 3-D reconstruction attempts require perfectly helical assemblies, which presents a serious limitation and confines the attempts to 15- or 16-protofilament microtubules, microtubule configurations that are very rare in nature. The rise of cryo-electron tomography within the last few years has now opened a new avenue towards solving 3-D structures of microtubule-MAP complexes that do not form helical assemblies, most importantly for the subject here, all microtubules that exhibit a lattice seam. In addition, not all motor domains or MAPs decorate the microtubule surface regularly enough to match the underlying microtubule lattice, or they adopt conformations that deviate from helical symmetry. Here we demonstrate the power and limitation of cryo-electron tomography using two kinesin motor domains, the monomeric Eg5 motor domain, and the heterodimeric Kar3Vik1 motor. We show here that tomography does not exclude the possibility of post-tomographic averaging when identical sub-volumes can be extracted from tomograms and in both cases we were able to reconstruct 3-D maps of conformations that are not possible to obtain using helical or other averaging-based methods.

A computational framework for ultrastructural mapping of neural circuitry.

Circuitry mapping of metazoan neural systems is difficult because canonical neural regions (regions containing one or more copies of all components) are large, regional borders are uncertain, neuronal diversity is high, and potential network topologies so numerous that only anatomical ground truth can resolve them. Complete mapping of a specific network requires synaptic resolution, canonical region coverage, and robust neuronal classification. Though transmission electron microscopy (TEM) remains the optimal tool for network mapping, the process of building large serial section TEM (ssTEM) image volumes is rendered difficult by the need to precisely mosaic distorted image tiles and register distorted mosaics. Moreover, most molecular neuronal class markers are poorly compatible with optimal TEM imaging. Our objective was to build a complete framework for ultrastructural circuitry mapping. This framework combines strong TEM-compliant small molecule profiling with automated image tile mosaicking, automated slice-to-slice image registration, and gigabyte-scale image browsing for volume annotation. Specifically we show how ultrathin molecular profiling datasets and their resultant classification maps can be embedded into ssTEM datasets and how scripted acquisition tools (SerialEM), mosaicking and registration (ir-tools), and large slice viewers (MosaicBuilder, Viking) can be used to manage terabyte-scale volumes. These methods enable large-scale connectivity analyses of new and legacy data. In well-posed tasks (e.g., complete network mapping in retina), terabyte-scale image volumes that previously would require decades of assembly can now be completed in months. Perhaps more importantly, the fusion of molecular profiling, image acquisition by SerialEM, ir-tools volume assembly, and data viewers/annotators also allow ssTEM to be used as a prospective tool for discovery in nonneural systems and a practical screening methodology for neurogenetics. Finally, this framework provides a mechanism for parallelization of ssTEM imaging, volume assembly, and data analysis across an international user base, enhancing the productivity of a large cohort of electron microscopists.

CTF determination and correction for low dose tomographic tilt series.

The resolution of cryo-electron tomography can be limited by the first zero of the microscope's contrast transfer function (CTF). To achieve higher resolution, it is critical to determine the CTF and correct its phase inversions. However, the extremely low signal-to-noise ratio (SNR) and the defocus gradient in the projections of tilted specimens make this process challenging. Two programs, CTFPLOTTER and CTFPHASEFLIP, have been developed to address these issues. CTFPLOTTER obtains a 1D power spectrum by periodogram averaging and rotational averaging and it estimates the noise background with a novel approach, which uses images taken with no specimen. The background-subtracted 1D power spectra from image regions at different defocus values are then shifted to align their first zeros and averaged together. This averaging improves the SNR sufficiently that it becomes possible to determine the defocus for subsets of the tilt series rather than just the entire series. CTFPHASEFLIP corrects images line-by-line by inverting phases appropriately in thin strips of the image at nearly constant defocus. CTF correction by these methods is shown to improve the resolution of aligned, averaged particles extracted from tomograms. However, some restoration of Fourier amplitudes at high frequencies is important for seeing the benefits from CTF correction.

Fibrils connect microtubule tips with kinetochores: a mechanism to couple tubulin dynamics to chromosome motion.

Kinetochores of mitotic chromosomes are coupled to spindle microtubules in ways that allow the energy from tubulin dynamics to drive chromosome motion. Most kinetochore-associated microtubule ends display curving "protofilaments," strands of tubulin dimers that bend away from the microtubule axis. Both a kinetochore "plate" and an encircling, ring-shaped protein complex have been proposed to link protofilament bending to poleward chromosome motion. Here we show by electron tomography that slender fibrils connect curved protofilaments directly to the inner kinetochore. Fibril-protofilament associations correlate with a local straightening of the flared protofilaments. Theoretical analysis reveals that protofilament-fibril connections would be efficient couplers for chromosome motion, and experimental work on two very different kinetochore components suggests that filamentous proteins can couple shortening microtubules to cargo movements. These analyses define a ring-independent mechanism for harnessing microtubule dynamics directly to chromosome movement.

Organization of interphase microtubules in fission yeast analyzed by electron tomography.

Polarized cells, such as neuronal, epithelial, and fungal cells, all display a specialized organization of their microtubules (MTs). The interphase MT cytoskeleton of the rod-shaped fission yeast, Schizosaccharomyces pombe, has been extensively described by fluorescence microscopy. Here, we describe a large-scale, electron tomography investigation of S. pombe, including a 3D reconstruction of a complete eukaryotic cell volume at sufficient resolution to show both how many MTs there are in a bundle and their detailed architecture. Most cytoplasmic MTs are open at one end and capped at the other, providing evidence about their polarity. Electron-dense bridges between the MTs themselves and between MTs and the nuclear envelope were frequently observed. Finally, we have investigated structure/function relationships between MTs and both mitochondria and vesicles. Our analysis shows that electron tomography of well-preserved cells is ideally suited for describing fine ultrastructural details that were not visible with previous techniques.

Automated electron microscope tomography using robust prediction of specimen movements.

A new method was developed to acquire images automatically at a series of specimen tilts, as required for tomographic reconstruction. The method uses changes in specimen position at previous tilt angles to predict the position at the current tilt angle. Actual measurement of the position or focus is skipped if the statistical error of the prediction is low enough. This method allows a tilt series to be acquired rapidly when conditions are good but falls back toward the traditional approach of taking focusing and tracking images when necessary. The method has been implemented in a program, SerialEM, that provides an efficient environment for data acquisition. This program includes control of an energy filter as well as a low-dose imaging mode, in which tracking and focusing occur away from the area of interest. The program can automatically acquire a montage of overlapping frames, allowing tomography of areas larger than the field of the CCD camera. It also includes tools for navigating between specimen positions and finding regions of interest.

The three-dimensional arrangement of intermediate filaments in Romney wool cortical cells.

The three-dimensional orientation and arrangement of intermediate filaments in Romney wool ortho-, meso-, and paracortical cells has been revealed using single axis high voltage electron tomography. Modelled tomograms confirm that intermediate filaments in orthocortical cells are arranged helically, with the helical angle progressively increasing from the centre to the periphery of macrofibrils. Intermediate filaments in meso- and paracortical cells display parallel arrangements differing mainly in packing density, with the mesocortex packed more tightly than the paracortex. The intermediate filament arrangements observed confirm expectations based on earlier two-dimensional transmission electron microscopy observations by the authors and other researchers. It is expected that these findings will contribute to a better understanding of the biological and structural basis of wool fibre curvature.

Three-dimensional ultrastructure of Saccharomyces cerevisiae meiotic spindles.

Meiotic chromosome segregation leads to the production of haploid germ cells. During meiosis I (MI), the paired homologous chromosomes are separated. Meiosis II (MII) segregation leads to the separation of paired sister chromatids. In the budding yeast Saccharomyces cerevisiae, both of these divisions take place in a single nucleus, giving rise to the four-spored ascus. We have modeled the microtubules in 20 MI and 15 MII spindles by using reconstruction from electron micrographs of serially sectioned meiotic cells. Meiotic spindles contain more microtubules than their mitotic counterparts, with the highest number in MI spindles. It is possible to differentiate between MI versus MII spindles based on microtubule numbers and organization. Similar to mitotic spindles, kinetochores in either MI or MII are attached by a single microtubule. The models indicate that the kinetochores of paired homologous chromosomes in MI or sister chromatids in MII are separated at metaphase, similar to mitotic cells. Examination of both MI and MII spindles reveals that anaphase A likely occurs in addition to anaphase B and that these movements are concurrent. This analysis offers a structural basis for considering meiotic segregation in yeast and for the analysis of mutants defective in this process.